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Abstract 

We study low-energy consequences of super symmetric 50(10) models with Yukawa unifi- 
cation ht = and /i;, = hr- We find that it is difficult to reproduce the observed mii/nir 
ratio when the third-generation right-handed neutrino is at an intermediate scale, especially 
for small tan [3. We obtain a conservative lower bound on the mass of the right-handed neu- 
trino Mjv > 6 X 10^^ GeV for tan /? < 10. This bound translates into an upper bound on 
the T-neutrino mass, and therefore on its contribution to the hot dark matter density of the 
present universe, U.^h?' < 0.004. Our analysis is based on the full two-loop renormalization 
group equations with one-loop threshold effects. However, we also point out that physics 
above the GUT-scale could modify the Yukawa unification condition hb = hr for tan /3 ^ 10. 
This might affect the prediction of mf,/m^ and the constraint on M^r. 
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1. Grand unification has been one of tlie leading guiding principles to build models of 
fundamental forces in nature. The introduction of a simple gauge group not only allows an elegant 
description of the electromagnetic, weak and strong forces, but also gives powerful insight into 
the structure of Yukawa couplings [@]. Supersymmetry, on the other hand, has provided a natural 
framework to solve the gauge hierarchy problem. On top of these theoretical ideas, there also 
appeared an experimental indication supporting supersymmetry combined with grand unification. 
Experiments as LEP |^ and SLC [§] have provided a very precise determination of the weak mixing 
angle. The resulting value has ruled out the minimal grand unified theories (GUTs), while showing 
a remarkable agreement with their supersymmetric versions based on the minimal particle content 
below the GUT scale Mq 0, 0- (The absence of particle thresholds between the weak and GUT 
scales is often referred to as "grand desert assumption".) In the same class of theories, the partial 
unification of Yukawa interactions {hf, = h^.) has also proven successful [0) i|) 01) while in the non- 
SUSY models it has not. These facts have strengthened the motivation for studying more detailed 
consequences of SUSY-GUTs with the grand desert. Recent work includes more theoretical efforts 
on GUT model building as well as more phenomenological studies such as those on cosmic neutralino 
abundance, proton decay and collider signatures. 

On the other hand, there are indications from astrophysics and cosmology that an intermedi- 
ate scale might exist in the middle of the proposed grand desert. One indication comes from the 
observed deficit in the solar neutrino flux [jlO|, which favors neutrino oscillations a la Mikheev- 



Smirnov-Wolfenstein (MSW) |Tl|] with non- vanishing neutrino masses of (9(10 ^) eV [|T2|. If we 



regard the neutrino oscillation as being of Ve^i^^i type, then the tiny masses required for the MSW 
effect suggest the existence of right-handed neutrinos having an intermediate mass of 10^°-10^^ GeV 
and acting as the source of a seesaw mechanism [|IH].r] Another indication comes from the den- 



sity fluctuations at large scales observed by COBE |14|, which have set the normalization of the 



Harrison-Zeldovich spectrum. The result, if one assumes that dark matter is only of the cold type, 



is in a moderate contradiction with the observations at smaller scales [p!5| , |T6|. It has been suggested 
that the existence of both cold and hot dark matter components solves this problem Jl^ . The only 
particle physics candidate for hot dark matter is a light neutrino with a mass around 10 eV. If we 
regard this neutrino as z/,-, then, again, a right-handed neutrino at an intermediate scale ~ 10^^ GeV 
is favoured. Moreover, it was also shown that scalar components of the right-handed neutrinos can 



play interesting cosmological roles in baryogenesis [T^ and inflation |[T9[| . 

Therefore, it is important to study the phenomenology of SUSY GUTs in which right-handed 
neutrinos are the only thresholds in the grand desert. To our knowledge, few studies of the low- 



energy phenomenology have been done within this context (see, e.g., H^])- Obviously, gauge 
coupling constant unification is not significantly affected because the right-handed neutrinos are 
gauge singlets. 

This letter studies the unification of the Yukawa couplings of the third family in SUSY GUTs 
with right-handed neutrinos at an intermediate scale. We base our analysis on 5*0(10) models with 

*We are assuming S'O(10)-like mass relations, where the neutrino Dirac masses are of the order of the up-type 
quark masses. 



1 



Yukawa unification ht = and hb = h^. We find that an intermediate scale right-handed neutrino 
has a critical influence on the predicted value of 171^/111^-. Consistency with the experimental value of 
mb/rrir implies a lower bound on the mass of the right-handed neutrino. This constraint translates 
into an upper bound on the r neutrino contribution fij, to the hot dark matter density in the 
present universe. Cosmologically interesting values of fl^, are strongly disfavored for small values 
of tan/3.[] 

2. We first describe the theoretical framework of our study. We assume a simple 5*0(10) GUT 
scenario where only small irreducible representations (irreps) like 10, 16, 16, 45, 54 and possibly 
singlets are present,[]and where matter fields come only in 16's. As for the Yukawa interactions, the 
large value of the top Yukawa coupling suggests that it originates from a renormalizable coupling. 
With our assumptions this can only be /igISsIOhIBs. (Even a 120 would not couple because of 
antisymmetry.) In general, the two light doublets Hu^d of the Minimal Supersymmetric Standard 
Model (MSSM) contain only fractions of those in IOh, due to mixing with doublets in other irreps. 
In our case these can be lO^'s, 16h's and 16h's, the latter ones being necessarily present in order 
to reduce the rank of the gauge group. The mixing strengths are in general different for and 
because SU{2)ji is broken at Mq- Then the GUT boundary condition on the Yukawa couplings is 
ht = h]\f = Suhc, hb = h^- = Sdhc, where Su,d are the Higgs mixing angles, and in general Su 7^ Sd- 
Here the meaning of ht, /iat, hb, and h^- is the obvious one. The particular has been 

studied before p5|, Actually the boundary condition 



ht = hN, hb = hr (I) 



is just a result of the underlying Pati-Salam SU{A) symmetry |2^, so that one may imagine sce- 
narios, maybe in strings, where eq. (P holds even without a unified simple gauge group. There are 
also more exotic scenarios with mixings between 16 and 10 matter fields |27|, while eq. (0) keeps 



holding. Moreover, there is a class of models where hb and hr originate solely from higher-dimension 
operators, while predicting the same relation We thus conclude that there is a wide class of 
models (5*0(10) and beyond) which lead to the boundary condition eq. (|l|). 

3. In order to illustrate the low-energy consequences of 6-r Yukawa unification when a right- 
handed neutrino is present, we first consider a simplified picture based on one- loop renormalization 
group equations (RGE), without threshold corrections. This also allows for a clearer comparison 

^ There is a claim po|] that 6-t Yukawa unification combined with the value of the top quark mass recently 
suggested by CDF excludes the region 1.5 < tan/? < 30. This claim, however, depends crucially on as{mz) — 0.125 
and mt ^ 174 GeV, for which the experimental support is not yet established. Thus, we take tan /3 as a free 
parameter in this letter. 

■i-Notice that large irreps drive the theory out of the perturbative regime; for example, if a 126 is present, the 
5*0(10) gauge coupling constant blows up below < 8Mg even with the minimum particle content 45 + 54 + 126 + 
126 + 3 X 16. It is also interesting to recall that, in superstring models, the above mentioned small irreps appear 
already at Kac-Moody level k = 2, whereas 126 can appear only for 5 < fc < 7 |Q. Phenomenologically, the 
presence of a 126 may also give an unacceptable relation hr — 3/ib [^. Finally, it is noteworthy that the absence 
of 126 and R-odd gauge singlets forbids right-handed neutrino masses at the renormalizable level, thus providing a 
dimensional argument for their being considerably smaller than Mq- 
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with the pure MSSM case, i.e. the case without right-handed neutrinos below the GUT scale. 
The RGE for the individual gauge and Yukawa couplings can be derived from the general formulae 
in Ref. both at one- loop order (needed here) and at two-loop order (to be used below), 

and we do not write them down explicitly. Here we rather focus on the evolution of the ratio 
R{fi) = hb{fJ^) / hrifJ^) from the GUT-scale value R{Mg) = 1 to the low-energy value R{mz)- The 
latter should be compared with the experimental value of mb/mT- scaled up to m^, which we denote 



by R^^^imz). We first recall the situation in the pure MSSM case |3^. The RGE for R reads 



cLR 
d^ 



R 
16^ 



hi + ?>{hl - h 



16 



-^3 



(2) 



where t = log/x. 
expression 



If all the Yukawa coupling constants are neglected, one obtains for R{mz) the 
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which is numerically larger than R'^^P{mz) by 30-80 %. The Yukawa contributions in eq. (0) 
are essential in reducing R{mz) from its 'pure gauge' value Rg to the observed value R^^P{mz). 
Actually, one could formally obtain arbitrarily small values for R{mz) by taking large enough ht 
at the GUT-scale [^. In practice, initial values of C(l) are required to reproduce the correct 
R'^^P[mz). On the other hand, if a right-handed neutrino with mass < Mq is present, the 
MSSM eq.(l) applies only for < /i < Mn- For Mn < < Mq the theory is the MSSM+N (in 
an obvious notation) and the RGE has a slight (but crucial) modification: 
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(4) 



The essential difference with the pure MSSM case is a manifest cancellation of all the Yukawa 
contributions, at least close to the GUT-scale, due to the boundary conditions ht = h^ and h^ = hr- 
As a consequence, even if one formally allows for arbitrarily large values of the Yukawa couplings 
at the GUT-scale, the low-energy value R{mz) cannot be made arbitrarily small, in contrast to the 
pure MSSM case. In particular, it may be impossible to reach the experimental value i?^^P(m^). 

To make this point more explicit, it is useful to integrate formally eqs.(^ and (^, and write 
R{mz) as 



R{mz) = Rg - e 



-Y 



where 



Y 



167r2 



hi{t)dt + / {hi{t) - hi,{t))dt + 3 / ihiit) - hi{t))dt 



(5) 



(6) 



Hz Jiff Jtz 

with Ia = logMyi. Notice that Y should be considered as a function of the initial values of the 
Yukawa couplings and the mass of the right-handed neutrino. Since all the three terms in Y turn 
out to be positive, R{mz) is smaller than Rg also in the present case. The question is whether or not 
Y can be large enough to reproduce the correct R'^^P{mz). Indeed the cancellation effect mentioned 
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above limits the size of Y , even for large initial Yukawa couplings. For a given < Mg, Y has 
a (finite) upper bound Ymax{.M]^), formally obtained by taking infinite Yukawa couplings at Mq- 
Since Ymax{,M]si) becomes smaller for smaller values of Mat, it may be impossible to get R{mz) as 
low as R^^^{mz) when Mj^ is lighter than a certain critical value. This is why a lower bound on 
Mjv appears. Actually, Ymax{,M]s[) crucially depends on tan/?. For fixed Mat, Ymax{,Mi^) is smaller 
when tan (3 is small, because and hj. are negligible and the last term in Y is missing. In fact this 
is the situation where stronger constraints on Mjq will be found. On the contrary, when tan/? is 
very large, the last term in Y plays a significant role (notice also the factor 3). Then, even if Mjv 
is light, Y can be made large enough as to reduce the value of R appropriately. 

The above qualitative discussion is supported not only by a numerical analysis, but also by an 
analytical argument which we sketch below. The argument establishes the existence of an upper 
bound on y as a function of M^r independent of the initial conditions on the Yukawa couplings. 
Even if one cannot solve exactly the coupled system of RGE, exact inequalities can be found 
that lead to the final result. We consider here the most interesting case where hb and are 
negligible. Then we can write Y = Yi + Y2, where Yi and Y2 are the first and second terms 
in eq. (^. We denote by hQ(t) the solution of the RGE for ht(t) in the pure MSSM case with 
boundary condition /io(Mg) — > 00. Its exphcit expression is hl(t) = 47r^(— F'(t)/3F(t)), where 

, 16 q 13 

F{t) = //«dt'(a3(MG)/a3(t'))'~ (a2(MG)/a2(t')) (ai(MG)/ai(t'))*- A study on the RGE of 
ht{t) for t < Mtv and t > Mn, with a given boundary condition htlMc), leads easily to the absolute 
upper bound ht{t) < hQ{t), valid for any t and independent of htiMc)- This gives immediately an 
upper bound on Yi, 

y,<yr(M.)4.og£i;M p) 

where we have replaced for convenience the arguments with Ma- A simple bound on Y2 can be 
derived e.g. by studying the RGE for the ratio R{t) = htif) /h]si{t), which reads 



1 dR 



R dt 167r2 



m - hi) - (y^3^ + ^/,) 



(8) 



If one formally integrates the above equation term by term between and tc? the right hand side 
contains Y2 and easily integrable gauge terms. The left hand side gives just — log-R(tAr), which is 
negative (one can easily prove that R{t) > 1 from a different form of eq.(^). Therefore one obtains 
an absolute bound 

Notice that the bounds (0) and (|) are independent of the initial values ht^Mc) = hN^Mc), so 
they hold even in the formal limit ht{MG) — 00. In conclusion, for a given Mjy < Mq we obtain 
explicitly an absolute upper bound Y^^^{M]s[) on Y which is a monotonically increasing function 
of M;v0 

Y < YmaxjMM) < Y:^,,{Mm) = Y*{Mm) + Y;{Mn). (10) 

§ Although YmaxiM^) is always monotonically increasing, the same applies to the analytic bound Y^^^{Mn) for 
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This in turn gives an absolute lower bound on R{mz) 

m,n,) > R, . e-'~<«») = R, ■ ( f^V* (^^f , (H) 
^ ' " \F(Mn)J \a3(MojJ \a,(Mo)J ^ ' 

proving that one cannot make R{mz) arbitrarily small even in the formal limit htlMc) — >■ oo. This 
bound can be compared with R^^P[mz)i so that a lower bound on can be inferred. 

4. We have shown an analytic bound on R in which /i^r was allowed to go to infinity at Ma- 
in practice /i^v will not be taken larger than 0{1). Then the effect of a right handed neutrino 
with mass around 10^^ GeV can be roughly estimated by treating it as a threshold correction to 
R at Mg 0. This gives R 1 + ^ logiMG/M^), which amounts to an 0(10%) increase in R. 



This 10% is a critical amount when comparing to the experimental value, since already in the pure 
MSSM case the predicted R lies in the upper part of the allowed range. Nonetheless it is clear that 
all other effects of the same order may have a huge impact on the bound on M^y. This is because 
the correction to R from the right-handed neutrino depends roughly logarithmically on M^. Then, 
in order to make our analysis more accurate, we have to discuss SUSY and GUT thresholds and 
use two-loop ROE. In fact the effects of these, and especially the first two, can possibly pile up to 
0(10%). We also attempt an estimate of the possible effects from physics above Mq, even though 
we are aware that our understanding of these can only be qualitative. Finally, by varying the QCD 
coupling as(m^) within its experimental range we get an 0(10-20%) effect, which we treat with 
the proper attention. Our reference value of asirnz) is the one from the Z-boson hadronic width 
otsiiT^z) = 0.124 ± 0.007 |3ll], and we allow a 2 cr range. 

To obtain quantitative constraints on the mass of the right-handed neutrino, we perform a 
numerical analysis based on two-loop RGE with one-loop threshold corrections. In order to be able 



to discuss the dependence on as{mz) P, 0, we take the following procedure. We define Mq as the 
scale where ai and a2 meet, and take ht = h^, hh = hr at the same scale.[] On the other hand, we 
take as{mz) as an independent parameter, and allow as ai = at Mq. This difference could 
be accounted for by GUT threshold effects. The gauge coupling constants in the MSSM at are 
a{^(mz) = 59.1 + f log ^2siisr + 1^ log a^Hruz) = 29.4 + f (4 log ^^^^^^^ + i log ^), and 

a^^{mz) = a~^{mz) + ^ log "^f^^^ ■ Here the logarithms take care of the matching between the 
SM and the MSSM; we take the additional Higgs doublet at and the rest of the SUSY particles 
at msusY- (The non- logarithmic parts and the translation from MS to DR are numerically small 
and have been neglected.) We also match the couplings in the MSSM-I-N to those in the MSSM at 

/i = MnE 

^fe(^Af)lMSSM=^fe(^Af)lMSSM+N ' ^r(MAr) |j^gg]y[= /l^(MAr) |^^gg^^_^_N (1 - £:), ^-^2) 

^t(^Af)lMSSM=^i(^Af)lMSSM+N " (^A^) ImSSM= (^A^) ImSSM+N ~ ^e). 



Mn ^ 10^° GeV only. However, this is just the region of interest. We add that stronger analytic bounds than the 
simple one described here can be found. 

^Numerically, Mg turns out to be around 2 • 10^'^ GeV. GUT threshold corrections on the Yukawa unification 
conditions will be discussed later. 

I'We define the right-handed neutrino mass Mjv by — MAr(Mjv), where Mn{^) is the DR running mass. 
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where e = ^^^2^ . h^i^j) in the MSSM is defined in such a way that the coefficient of the 
dimension-five operator {LHuY is h%{fi)/MN- 

For any fixed values of (/if, hh) at Mq, the actual calculation is done as follows. We use the 
two-loop RGE of the MSSM and the MSSM-I-N. By using an iterative procedure, we find the values 
of Mg, ai^Mc) and asi^Mc) which reproduce the values of ai(m^), a2{mz) and as{mz) defined 
above. At the same time we obtain the values of the Yukawa couplings at mz- The resulting 6-r 
Yukawa ratio -R("^z)Imssm then matched to the ratio of the corresponding masses in the broken 
electroweak theory R{mz) = mi,{mz) /rriT-imz), 

R{mz) = {I + h - K + fn) R{mz)\MssM > (13) 

where kf,, k^- are the threshold corrections due to the SUSY particles, while fn is that from the 
additional Higgs doublet in the MSSM. The exact expressions are found in Ref. [^. In the 
following we will only focus on the logarithmic terms in /c^,-. We will discuss the potentially 
large non-logarithmic ones (fc^ in Ref. [^) at the end of the letter. Note that also -R("^z)Imssm 
depends implicitly on msusY and due to the threshold dependence of the MSSM gauge coupling 
constants a^^{mz) (see above). For instance, in the small tan/3 region, the dependence of R on 
the SUSY particle masses turns out to be approximately^ 

^ - 7^ -0-1 log + 0.1 log . 14 

R ZTT V mz mz) 

The threshold corrections roughly cancel when vnsusv = which is consistent with the results in 
Ref. p. In the final step, the predicted R{mz) is compared to R^^'^ijnz)- Notice that also R^^^imz) 
has an important dependence on ^^(m^), since it is obtained by scaling the 6-r mass ratio up to 
niz- GUT-scale threshold corrections could also affect the boundary condition in eq. (|1]). Though 
we neglect such corrections in the following numerical analysis, we will discuss their possible effects 
later. 

We first illustrate the dependence of R on ht{Mc) for the region tan/? < 10, where both hb-, 
hj. get renormalized homogeneously and the result depends very weakly on tan /?. We show curves 
for different values of M^r in Fig. 1, taking msusY = m^ = mz, = h^- = 0.01 at M^, and 
as{mz) = 0.11 (Fig. la) or 0.12 (Fig. lb). For comparison, we also show the values of R^^^{mz) 
corresponding to the MS mass m;,(mb) = 3.9, 4.15, and 4.4 GeV. This is the range which is obtained 
from QCD sum rules |3^, ^.[^ First of all, it is clear from the figures that it becomes harder to 
reconcile R{mz) with R^^P{mz) for larger as{mz) and lower M^. For ^^(m^) = 0.12, one needs 

** In a more general case where the SUSY spectrum is nondegenerate, the dependence of R on the SUSY particle 
masses can again be summarized by an effective msusv ■ It turns out that the dependence of R on colorless particle 
masses is very weak, and msusv can be interpreted as a geometric average of squark and gluino masses. After 
considering generic splittings in the SUSY spectrum, we found that the value of R varied only between two extreme 
cases, msusv = 1 TeV, mA = mz, which gives the smallest R, and msusv = mz, mA = 1 TeV which gives the 
largest R. We stress that the approximate formula eq. ( p^ cannot be applied for mA ^ mz- 

^^The uncertainty on mi,{mi,) is dominated by our lack of knowledge on the 0{a^) corrections to QCD sum rules. 
See Ref. il. 
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very large htiMc) ^ 3 even for Mn = Mc-f^ Moreover, even if increasing ht{Mc) makes R{mz) 
decrease, the suppression effect becomes weaker as we lower M^v- Note also that the value of 
R{mz) depends only weakly on ht{Mc) for ht^Mc) ~ 2, especially for < 10^"^ GeV. Such 
a 'fixed point' behaviour and the dependence on Mat were expected on the basis of the analytic 
discussion above (see e.g. eqs. (p!0| - |Tl]) ). As we lower Mat, we can check whether the curves reach 
the region of R^^'^{mz) and then infer a lower bound on M^r. In the following, when scanning the 
space of (ht^Mc), hh^Mc)), we will only take for definiteness htfi^Mc) < 2. This reference value is 
motivated both by the above observation and by perturbativity reasons (see also point 2 below). 

Next we show lower bounds on Mjv as functions of tan f3 in Fig. 2. In order to obtain conservative 
bounds we take the maximum value of mb{mh) = 4.4 GeV. For most of the tan/3 values, we obtain 
stringent lower bounds on M^. We show the cases of two representative SUSY particle spectra 

(a) msusY = 1 TeV (conservative) and {h) nisusY = (more stringent), while always keeping 
niA = TTT-z (conservative). Curves are shown for three values of ^^(m^) = 0.110 (solid), 0.117 
(dotted), and 0.125 (dashed), where the lower ones correspond to case (a) and the upper ones to 

(b) . The values of M^r above the curves are consistent with 6-r Yukawa unification. The constraint 
becomes stronger for larger values of as{mz), smaller values of msusY, and larger values of m^. The 
curves do not extend to the region tan/5 ^ 58 because of the constraint = h^- < 2 at Mq- When 
the curves reach M^, it means that the lower values of tan/5 are not consistent with 6-r Yukawa 
unification even when = Mq (with htiMc) < 2). Possible effects from the non- logarithmic 
SUSY threshold corrections and physics above the GUT-scale will be discussed in points (1) and 
(2) below. 

The lower bound on M^v has a very interesting implication for the value of fl^, the u^- contribution 
to the hot dark matter density of the present universe. The mass of u^- is related to Mj^ and hj^ 
via the seesaw formula,^] 

Mj, ' ^^^^ 

where v = 174 GeV and /iat in the effective theory below Mn was defined after eq.(|12]). On the other 
hand, the cosmic energy density of a light neutrino of mass rrii, ^ 1 MeV is simply proportional to 
m 



33 



Q^h' ^ ^'^ . (16) 
91.5 eV ^ ^ 

Here, h is the normalized Hubble constant, h = /fo/(100 kmsec"^ Mpc~^). Note that hN{fnz) 
cannot be larger than an infrared fixed point value, which is h]\]{mz) < 0.8 for the cosmologically 
interesting region M^v ^ 10^^ GeV. Using that value in eqs. (^^-|T6|), we obtain the maximum 
possible value of VLj^h"^ for a given M^y. By requiring that Vj. gives a certain contribution to f2, one 

MPor as{mz) ^ 0.12 and htiMc) < 2, the predicted value of R{niz) is larger than the experimental upper bound. 
Consistency could be possibly restored by allowing 0(5%) GUT- or SUSY-scale threshold corrections. See discussion 
below eq. (19). 

*Here we neglect SUSY- and weak-scale threshold corrections on m,^^, which are at most 0(5%). Neglecting 
them is also justified because they affect the constraint on Mj^ only linearly. In contrast, the analogous corrections 
to R{mz) discussed above are larger and affect the constraint on A/jv exponentially. 
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can put an upper hound on M^. We indicated such bounds for Vtyh'^ = 0.01, 0.1, and 1 in Fig. 2 
with dotted hues. We find that the cosmo logically interesting range Q^h"^ ^ 0.1 survives only for 
large tan/?, especially for larger values of as{mz). For as{mz) ~ 0.125, we find solutions only for 
very large tan/5, and they completely disappear for as{mz) > 0.13. In other words, this analysis 
puts an upper bound on the cosmic hot dark matter density as a function of tan/?. For instance, 
we obtain 



for tan/? < 10 even with the most conservative parameters as{mz) = 0.11 and case (a). 

5. We now discuss the GUT-scale threshold corrections on ht = h^, ht = hr- They mainly 
come from the mass splittings within the heavy gauge multiplet as far as there are no other sizable 
Yukawa couplings of I63 to other GUT-scale fields ||2^. For instance, mass splittings of (9(10) 
modify the boundary condition by less than 5 A 5 % correction on ht/hi\f leads to less than 
1 % change in the prediction of R{mz) so that it will not change the constraint at all. However 
a 5 % threshold correction on /i(,//ir results in a 5 % correction in R{mz), leading to a two order 
of magnitude change on Mjv (see Fig. 1) and one might imagine a situation where the combined 
effects of as{mz) = 0.11, SUSY and GUT thresholds pile up to allow = 10^^ GeV, even at small 
tan/?. Situations like this are however rather extreme. When as{mz) — 0.12 (close to the centre of 
its range) and SUSY particles are light (as they should be), then GUT threshold corrections must 
amount to 15% to make = 10^^ GeV possible, which requires (9(1000) mass splittings in the 
gauge multiplet. Relaxation of the maximum allowed htlMc) up to 3.3 could also weaken the 
constraint on VL^h'^, but no more than by a factor of 7. 

Our bound on VLy has been derived in the approximation in which family mixings are neglected. 
Then one might wonder whether, or under what conditions, the mixings can sizably affect the 
bound. In order to evade our bounds, the full mass matrices should satisfy the following two 
requirements: i) The interaction eigenstate C I63 with Yukawa coupling = ht has an 0{1) 
overlap with a mass eigenstate of mass ~ Mq- Then the effective appearing in the RGE below 
Mg is smaller than ht, and the effect of ht in the running of R is no longer "exactly" compensated, ii) 
There still is a left-handed neutrino which contributes to a sizable portion of f2. In order to study the 
consequences of these two requirements, we consider the case of two families, which may represent 
the second and third ones. The superpotential reads as = hl^LiCjIId + H'^LiNjHu + ^MijNiNj 
with i,j = 2,3, where L, e and are, respectively, lepton doublet, right-handed charged lepton, 
and right-handed neutrino superfields. The smallness of CKM angles and 5*0(10) relations among 
Yukawa matrices suggest that both h]^ and /i^ are hierarchical in the same basis. By going to the 
basis where /iat is diagonal, we parametrize 



^^As explained above, our procedure implicitly requires a GUT-scale threshold correction on 0:3 (Mfj) to reproduce 
as{mz)', however, the required correction amounts to at most 5 % if we vary as{mz) = 0.11-0.14 and msuSY = mz- 
1 TeV. Though the corrections on and hh/hr are not directly correlated, this value gives us a rough idea of the 
magnitude of the GUT-scale threshold corrections. 



n^h'^ < 0.004, 



(17) 




(18) 
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where Mij ^ (9(Mg) and we expect e ~ mc/mt ~ 10^^. Requirement i) implies that at least one 
of the Mij should be 0{Mg), but not M22 alone. The light neutrino mass matrix can be written as 



22 



Requirement ii) is expressed as (m^ v'^siv? (3/{5 ■ Mq), where mi^2 are the eigenvalues 

of A4i^, and we need 6 ^ 10~^ for a neutrino mass in the eV range. Therefore one needs 

I DetM I < 6 (e^M^g + 2e^M^r^ + M^^ ^'"^ Mg, (21) 

where we used /itv ~ 1- The above inequality implies a hierarchy in the eigenvalues of M. This leads 
to strong constraints on M. If all the entries of M are 0{Mg), we need a fine-tuning of 0{6) to 
obtain the small determinant in eq. (|2l|) f\ We can avoid such a fine-tuning only if M is hierarchical 
in the same basis where h is diagonal with M22 ~ Se'^Mc, M23 ^ S^^'^eMc and M33 ~ Mq- In this 
case M could be viewed as a simple consequence of an abelian horizontal symmetry, but notice that 
the expansion parameter S^^'^e would be rather small. Moreover, the hot dark matter candidate 
would be predominantly z/^, and the MSW oscillation should occur between z/g and Ur- This is a 
new possibility which may be worth of further study. Nonetheless both the latter possibility and 
the fine-tuned case mentioned above show rather extreme features, thus enhancing the importance 
of the cases to which our analysis correctly applies. These cases include the reasonable situation in 
which there is no hierarchy among right handed neutrino masses, and they all decouple from the 
theory essentially at a single scale. 

6. Finally we point out two generic uncertainties in 6-r Yukawa unification which might affect 
the predictivity on R. These exist even in the pure MSSM case. Of course, such uncertainties may 
also affect the bound on M^. 



[1) SUSY-scale non-logarithmic threshold corrections to R for large tan/3 p4|, |3^. Typically the 



largest corrections appear in the b mass via two diagrams, one involving gluino propagation 
and the other involving higgsinos. The resulting correction to ruf, can be written as 

^ = t,^p(?^^ + Ji.t^), (22) 

where rrig and At are respectively the gluino mass and the stop trilinear coupling, while mi 
and 7712 represent effective SUSY masses out of the loop integrals. In particular mi roughly 



*The tuning implies that that the leading contribution to M is rotated by 0(1), with respect to the leading one 
in /ie_jv- It may be interesting to ask whether this can result from a flavour symmetry. 
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corresponds to the maximum between the sbottom and gluino masses, while m2 represents 
that between stop and higgsino (more exact expressions are found in [0]). When the SUSY 
parameters are all of the same order of magnitude, the typical size of the expression inside 
brackets in eq. (|22|) is ~ 1%. Then for tan/5 = 0{1) we can safely neglect them. On the 
other hand for tan/5 ~ 10 their effect could be a 10% reduction of R (notice that the sign 



of eq. (p^) is not fixed). As already stated, a 10% reduction in R is critical for the bounds 
on Mat, so that ~ 0.1 could be allowed for tan/5 > 10. Notice however that 6mh depends 
rather strongly on the features of the SUSY spectrum and becomes negligible when either fj, 
or the gluino mass (and At) are somewhat smaller than the squark masses. Only after an 
experimental determination of the SUSY mass parameters, will we be able to make conclusive 
statements on the tan/? ^ 10 region. It remains however true that in a relevant region of 
parameter space (with small Snih), the bounds shown in the figure are valid.f] 

(2) Corrections from physics above the GUT-scale. In our analysis we allowed a relatively large 
htlMc) < 2, which implies that the Yukawa Landau pole, and presumably some new physics, 
can be as close as SM^.! In general, the presence of new physics at a scale A > Mq induces 
non-renormalizable operators in the effective theory at the GUT-scale. Such operators will 
in general affect eq. (|^) with contributions scaling as powers of M^/A = rj, and one could 
expect 7] ~ 0.1 {e.g. from the above observation). Indeed, the order of magnitude of this 
ratio naturally suggests that flavour mixings and mass ratios are generated by the same type 
of effects. From this viewpoint, it is legitimate to make the order of magnitude interpretation 
rj ^ Vcb — 0.05. Flavour scenarios of this type have been depicted in Ref. We 
will therefore use the indicative value above to estimate the effects of non-renormalizable 
operators, which we classify below. 

A first class of operators is obtained by inserting 5*0(10) breaking Higgs fields in the renor- 
malizable term I63IOHI63. Our assumption that there be small irreps only (45's, 54's, 
16's, 16's and lO's ) constrains the lowest possible correction to be 0{rf) ~ 1%. These 
might come from 16310h(45)(45)163/A2 and 16310h(16h) (16h)163/A2.0 Such effects can 
be safely neglected. 

A second class of operators involves the other Higgs fields 10^, 16h and 16h, rather than the 
IOh appearing in the renormalizable term. The only possible dimension five operators are 
16310^1(54)163, 16310'h(45)163, 16316h(16h)163 and 163T6h(T6h)163. None of them 



^ It must be added that the approximate symmetries which render large tan [3 natural suppress Snu, [P4[ . On the 
other hand, these approximate symmetries require rogugy ^ P'^'z because of the LEP constraint /i. mi/2 raz- 
Therefore the casejirrig / rngjj ~ 1 ^"^^ the case l^rnglrn^g^gY ~ 1/ tan/3 require a comparable fine-tuning in the 
Higgs sector |3^, However, phenomcnological constraints from the observed rate for 6 ^ 57 again seem to favor 
the small Srrib case [|6[ 0. 

tif one requires perturbativity up to Mpianckl ^ one needs htiMo) < 1.5. On the other hand, notice that the 
value ht{MG) = 3.3 allowed in previous studies |^ implies that the Landau pole is closer than 2Mg. 

§Note that (54) preserves 0(6) x 0(4) (Pati-Salam group), and cannot modify eq. (1). Also, 16310h(45)163 
respects eq. (1) due to the symmetry between two I63. 
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affects hb = hr, while the last one may induce a correction 0{ri) ~ 10% to ht = Hn- Such a 
correction does not affect our analysis, as we already stated (see discussion below eq. (17)). 
On the other hand, dimension six operators can affect ht = hr. Actually, the relevance of 
this class of operators strongly depends on the magnitude of hb, hr themselves, i.e. on tan/?. 
When hb^r ~ V (or tan/3 ^ rjmt/mb ~ 10), dimension six terms can give corrections which 
are Oirf/hb^r) ~ ~ 5-10%. We recall that, for hb^r ^ 1, sits mainly in other Higgs 
fields than IOh- In fact, for hb^r ~ Vi the above dimension five operators become a natural 
source for these Yukawa couplings. When hb^r ~ tan/? < 10, the dimension five 

operators themselves have to be suppressed, since they would typically yield 0{ri) Yukawa 
couplings. This could result from a flavor symmetry. In this situation the dimension six 
operators become relevant and can finally even dominate hb^r- Since these terms can give 
both hb = hr and 3hb = —hr (via "composite" 126 combinations), or a mixture of both, the 
only way not to lose predictivity completely is to suppose that only the first type of terms 
exists. With this assumption the corrections to hb = hr are again expected to be at most 
0(7]), and could be actually 0(77") (with n > 1) in more specific flavour models. 

In short, hb = hr is an automatic consequence of gauge symmetry and field content alone for 
hb,T ~ V- For smaller hb,r, the relation hb = hr should probably result from flavor symmetries. 

We are aware that the last class of corrections is less under our control than all the other 
ones discussed so far. In fact, an opposite viewpoint would emerge if neutrino-oscillation 
experiments like CHORUS or NOMAD should find m^^ ^ eV. This discovery could indeed be 
used to probe the existence of higher dimensional effects. For instance, should the observations 
suggest Mn = 10^2 GeV, then we would need R{Mg) = 1.05, 1.13, 1.23, for a,(mz) = 0.110, 
0.117, 0.125, respectively. It is manifest that for the larger values of as we need the second 
class of effects (and larger than expected), while for the smaller as GUT thresholds alone 
could account for the deviation of R{Mg) from unity. These remarks could be a useful guide 
in building realistic GUT models. 

Therefore, there are possible corrections to the prediction of R{mz) for both tan/? ^ 10 and 
tan/? ^ 10 from very different origins (indeed the reason is the same; namely nib ^ f^t)- For 
large tan/3, non-logarithmic SUSY threshold effects can modify the prediction, but they can be 
calculated after the SUSY spectrum will be known. For small tan/?, corrections arising from 
physics beyond the GUT-scale could lead to two different classes of effects. When only first class 
effects are present, like in Ref. [0, then the condition in eq. (Q) is robust and so are our bounds 
on Mtv. On the other hand, when tan/3 <^ mt/rrib, reasonable flavor physics scenarios with the 
second class of operators may generically lead to a (9(5-10%) shift in the boundary condition, 
which is typically non-negligible for our purposes. However, in specific models this effect may also 
be further suppressed by symmetries. Moreover, even with (9(5-10%) corrections, there still are 
stringent limits on for asimz) ^ 0.12. Note also the new physics scale A could be larger 
> lOOMc when ht^Mc) ^ 1.5. At large tan/?, GUT-scale uncertainties are under control, so that 
the measurement of SUSY parameters will be enough to make our analysis very accurate. 
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7. In summary, we have studied the impact on the low energy value of mb/mr of an intermediate 
mass right-handed neutrino in a class of SO{10) models. The analysis generally results in a lower 
bound on the mass Mn of the right handed neutrino, and thus in an upper bound on the r neutrino 
mass and on its contribution fij, to the hot dark matter density in the present universe. In order to 
do so we have performed a two-loop study of the RGE and discussed threshold effects at the GUT 
and SUSY scales. Depending on the values of and tan/3 the bound varies from a very strong 
one to a weaker one. 

i) When tan/3 ^ 10, there are uncertainties in the bound on M^v from the yet-unknown non- 
logarithmic SUSY threshold corrections. When these effects are maximal, then is allowed to 
span all the phenomenologically interesting region (down to Mn ~ 10^2 Qgv). On the other hand, 
when these effects are small, then only for very large tan/? > 50 can we reach the cosmologically 
interesting region Mn ^ 10^3 GeV. We stress that effects at the SUSY-scale threshold will be 
known once SUSY particles are discovered and their parameters measured, so that such effects will 
be no longer ambiguities. Instead, we will be able to evaluate the corrections, which might even 
make the constraint stronger. 

ii) For small tan/5, we get rather stringent bounds on Q^. For as{mz) = 0.11,0.117,0.125 we 
respectively get fl^h'^ < 4x 10"^, 3x 10~^, 6x 10"^ for tan j3 < 10. Corrections from SUSY thresholds 
do not modify these bounds. If we allow for "maximal" GUT threshold effects of (9(5%), these 
bounds can go up by about two orders of magnitude. We point out that there might be additional 
corrections to hb/hr = 1 at Mq from physics beyond Mq- Assuming that the new physics is 
responsible for the flavor structure, we estimate the size of the corrections to be 0{Vcb) ~ 5%. Even 
taking all these possible ambiguities into account, the cosmologically interesting region ^ 0.1 
is allowed typically only for as{mz) ^ 0.12. Improvement in the experimental knowledge of 
[p9| will allow tighter bounds, especially if the value will converge towards the present central one 
>0.12. 
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Note Added 



After completing the main part of the work described in this letter and presenting it at a conference 
|D|, we received a preprint by Vissani and Smirnov discussing a similar topic [ET|. Their basic 



conclusion is that an intermediate mass right-handed neutrino disfavors 6-r Yukawa unification, or 
viceversa, which is the same as ours. However when deriving quantitative lower bounds on Mtv, 
their analysis is based on one-loop RGE with tree-level matching, while ours is based on two-loop 
RGE with one-loop matching. We consider this to be necessary due to the high sensitivity of the 
bounds on to (9(5-10%) effects which might arise from SUSY-scale and GUT-scale threshold 
corrections. In addition, they do not discuss the possible relevance of points (1) and (2) we discussed 
above. However, in the cases where all such effects are negligible, our result is consistent with theirs. 
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Figure Captions 



Fig. 1 An illustration of the effect of the right-handed neutrino on R{mz) = mb{mz) /mr{mz) ■ The 
curves show the dependence of R{mz) on ht{MG) = hNiMc), for the values of — Mq, 
10^5 GeV, 10^^ GeV, lO^^ GeV, lO^^ Qgy, and 10^° GeV. The experimentally allowed values of 
R{mz) are also shown. Other input parameters are taken as msusY — f^A — f^Zi hb{MG) — 
hriMo) = 0.01, and as{mz) = 0.11 (a) or 0.12 (b). 

Fig. 2 Lower bounds on the right-handed neutrino mass Mjv from b-T Yukawa unification. Curves 
for three values of ttg = 0.110 (solid), 0.117 (dots) and 0.125 (dash) are shown, both for the 
SUSY particle spectra (a) msusY = 1 TcV, = mz (lower) and (b) msusY = = 
(upper); see text. The possible cosmic energy density ili.h'^ of the r-neutrino hot dark matter 
is also shown for comparison in dot-dashed hues. 
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